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ABSTRACT 

We present a measure of the hard (2-8 keV) X-ray luminosity function (XLF) of Active Galactic 
Nuclei up to z ~ 5. At high redshifts, the wide area coverage of the Chandra Multiwavength Project is 
crucial to detect rare and luminous {Lx > 10^^ erg s~^) AGN. The inclusion of samples from deeper 
published surveys, such as the Chandra Deep Fields, allows us to span the lower Lx range of the 
XLF. Our sample is selected from both the hard {z < 3; f2-8 kcv > 6.3 x 10~^^ erg cm~^ s~^) and 
soft {z > 3; /o. 5-2.0 kcv > 1-0 X 10~^^ erg cm~^ s~^) energy band detections. Within our optical 
magnitude limits (r', i' < 24), we achieve an adequate level of completeness (> 50%) regarding X-ray 
source identification (i.e., redshift). We find that the luminosity function is similar to that found in 
previous X-ray surveys up to 2; ~ 3 with an evolution dependent upon both luminosity and redshift. 
At z > 3, there is a significant decline in the numbers of AGN with an evolution rate similar to that 
found by studies of optically-selected QSOs. Based on our XLF, we assess the resolved fraction of the 
Cosmic X-ray Background, the cumulative mass density of Supermassive Black Holes (SMBHs), and 
the comparison of the mean accretion rate onto SMBHs and the star formation history of galaxies as 
a function of redshift. A coevolution scenario up to 2; 2 is plausible though at higher redshifts the 
accretion rate onto SMBHs drops more rapidly. Finally, we highlight the need for better statistics of 
high redshift AGN at z > 3, which is achievable with the upcoming Chandra surveys. 
Subject headings: galaxies: active — quasars: general — X-rays: galaxies — surveys 



1. INTRODUCTION 

Our present understanding of the evolution of accret- 
ing SMBHs over cosmic time comes from the measure 
of the luminosity function (i.e., the number undergoing 
a detectably luminous phase within a specific co-moving 
volume as a function of luminosity and redshift) of Active 
Galactic Nuclei (AGN) . Energy production through mass 
accretion onto SMBHs allows us to observationally iden- 
tify these sites as the familiar AGN with Quasi-stellar 
objects (QSOs) the most luminous example. Hence, the 
AGN luminosity function provides a key constraint to 
discern the underlying physical properties of the popu- 
lation (i.e., black hole mass and accretion rate distribu- 
tions as a function of redshift) and thereby elucidate the 
mechanisms (i.e., galaxy mergers and/or self- regulated 
growth) that are instrumental in their formation and evo- 
lution. 

To date, an enormous effort has been undertaken to 
measure the luminosity function over the wide range 
in luminosity spanned by AGN out to high redshift. 
The bright end has been well established to z ^ 5 by 
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optical survevs (e.g.. [ Richards et aT]l2006t ICroom et al.l 
[2OOI LWolf et al."2003) which primarily select QSOs us- 
ing a multi-color photometric criteria. The most dra- 
matic feature found is the rise and fall of the co-moving 
space density with peak activity at z ~ 2.5. With an 
unprecedented sample of over 20,000 QSOs in th e 2dF 
QSO Redshift Survey f2QZ). ICroom et al.l (jlOO?) con- 
vincingly show a systematic decrease in luminosity (pure 
luminosity evolution; PLE) from z = 2 to the present, 
in agreement with past surveys (e.g.. 



Schmidt fc GreenI 



119831 : iBovle et al.lll988t iHewett et aLlll993D . which find 
very few bright QSOs in the local universe. This fading 
of the luminous QSO population is thoug ht to be due to a 
decre ase in the mass accretion rate (e.g.. ICavaliere et"ari 
I2OOOI ) that appears to be intimately related to the 
order-of-magnitude decline of the c osmic star forma- 
tion rate from z ^ 1 to the present fBovle fc Ter levichI 
[1998; Franccsc hini et al.Ll999 : Merloni ct al. 2004) The 
dropoff in the space density at z > 3 fOsmen 11982 



Warren et al.lfT99l[Sctmridt et al.lfT995t iFan et al.ll2001 



Wolf et al.ll2003[ r may be indicative of either the detec 



tion of the onset of accretion onto young SMBHs or 
a high-redshift population that has b een missed, pos- 
sibly under a veil of obscuration (e .g., [Alexander et al.l 
l2005t[M"artinez-Sansigre et al.ll2005[ ). Excessive amounts 
of dust and gas may be ubiquitous in galaxies at 
early epochs due an increase in the merger rate 
(jKartaltepe et al . 2007 ) tha.t indu ces high star formation 



rates (e.g.. IChapman et al.ll2005l ) 



It has been evident for quite some time that optical 
surveys miss a significant fraction of the AGN popula- 
tion. They fail to find the majority of AGN due to a 
steeply declining luminosity function with the low lu- 
minosity end severely affected by host-galaxy dilution. 
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Though current optical selection techniques do show con- 
sider able improvement (|Richards et al.l[2005t iJiang et all 
|2006() . they still fail to account for many low luminosity 
AGN. Of equal significance, many AGN (e.g., Seyfert 2s, 
narrow line radio galaxies) are missed due to dust ob- 
scuration (causing their optical properties to differ from 
the type 1 QSOs) and can o nly be adequate ly selected 
in the low redshift universe (|Hao et al.ll2005[ ) based on 
their highly ionized, narrow emission lines. Fortunately, 
AGN can now be efficiently accounted for by selection 
techniques in other wave bands such as the X-ray as 
demonstrated in t his work, radio (e.g. , IWall et 31112005( 1 
and infrared (e.g.. lPolletta et aLll2006l ). As further elab- 
orated below, current models based on recent observa- 
tions continue to attribute the bulk of the Cosmic X-ray 
Background (CXRB), the previously unresolved X-ray 
emission, to these various types of obscured AGN. 

Over more than two decades. X-ray surveys have been 
improving and extending the known AGN luminosity 
function by including sources at low luminosity, with or 
without optical emission lines, and hidden by a dense 
obscurin g medium. The Extended Medium Sensitivity 
Survey (jGioia et al.l Il990( ) was one of the first surveys 
to measure the X-ray luminosity function (XLF there- 
after) out to the QSO peak using a sample of just over 
400 AGN detected by the EINSTEIN Observatory. 
Since the survey only probed the more luminous AGN 
[Lx > lO^'' erg s~^) above z = 0.3 due to the bright 
flux limit (/0.3-3.5 kcv > 5 x 10~^^ erg cm~^ s ~^), i t was 
quite u nderstandable that iMaccacaro et al.l ()1991l ) and 
Idella C oca et al. (1992) found the XLF to behave simi- 
larly to the optical luminosity function (i.e., PLE) with 
a decreasing space density from z ~ 2 to the present. 
The ROSAT satellite with its increase in flux sensitiv- 
ity (/0.5-2.0 kov > 10~^^ erg cm~^ s~^) enabled AGN 
to be detected at lower luminosities, and out to higher 
redshifts (z ^ 4.5). Surveys ranged from the wide area 
and shallow ROSAT Bright Survey (f^m ~ 10"^^ erg 
cm-2 s-i; 20000 deg^: ISchwope et al.ll20M ) to the deep 
and narr ow Lockman Hole (ium ~ 10~^^ erg cm~^ s~^; 
0.3 deg^; iLehmann et al.l [20011 ). With a compilation of 
690 AGN from these fields an d other available ROSAT 
survevs (iBower et al.' Il9 96t lAppen zeller et al.l 119981: 
McHardv et al. 1998; Zamo rani et al. , 1998: Mason et all 



2000( ) tha t effectively fiU in t he parameter space of flux 
and area, iMivaii et al.l (|2000f ) were able to resolve « 60- 
90% of the soft CXRB into point sources and generate 
a soft XLF that extended beyond the QSO peak. They 
found that the XLF departed from a simple PLE model, 
now well known to describe X-ray selected samples as fur- 
ther elaborated below. A luminosity-dependent density 
evolution (LDDE thereafter) model was required due to 
the slower evolution rate of lower luminosity AGN com- 
pared to that of the bright QSOs. The limited sky cover- 
age at the faintest X-ray fluxes achievable with ROSAT 
prevented an accurate measure of both the faint-end 
slope at z > 1 and the overall XLF at high redshift since 
only a handful of AGN were identified above a redshift 
of 3. 

In the current era of Chandra and XMM-Newton, X- 
ray surveys are now able to detect AGN and QSOs not 
only enshrouded by heavy obscuration but those at high 
redshift (z > 3) with greatly improved statistics due to 
the superb spatial resolution and sensitivity between 0.5 



to 10 keV of these observatories. Previous X-ray mis- 
sions such as EINSTEIN and ROSAT as described 
above were limited to the soft band which biases sam- 
ples against absorpti on. The ASCA observatory (e.g., 
'Akiva ma et al.l [2003) successfully found many nearby 
absorbed AGN but lacked the sensitivity to detect the 
fainter sources contributing most of the 2-8 keV CXRB. 
With deep observations of the Chandra Deep Field North 



jep 

( CDF-N; D 3randt ct al. 2001), Deep Fiel d South (CDF- 



S; Rosati ct al. 2002) and Lockman Hole fHas inger et al 
| 2001) . a large frac tion, between ~ 70% (Worslev et al 
I2OO5I ) and ~ 89% ()Moretti et al.ll2003[ ) of the hard (2-8 
keV) CXRB has been resolved into point sources. Many 
of the hard X-ray sources found so f ar arise in optically 
unrem arkable bright galaxies (e.g., iBarger et al.l [20031 : 
iTozzi et al. 20Q1|), which can contain heavily obscured 
AGN. 

A more robust luminosity-dependent evolutionary 
scheme has emerged from recent measures of the XLF. 
With the inclusion of absorbed AGN from Chandra and 
XMM-Newton surveys, lower luminosity AGN are clearly 
more prevalent at lower redshifts (z < 1) than those of 
high luminosity that peak at z ^ 2.5. This behavior 
is due to the flattening of the low luminosity slope at 
higher redshifts that has been well substantiated with 
hard (2-8 k eV) X-ray sel ected surveys (Cowie et al. 200j; 
Ueda et al. 2003; Fiore et all [20 03: Bargcr ct ali 120051 : 



Silverman et al... 2005b; iLa Franca et al..,2005i) . Using a 



highly complete soft (0.5-2.0 keV) band sel e cted sam- 
ple of over 1000 type 1 AGN, iHasinger et aH (|2005( ) has 
shown that this LDDE model accurately flts the data 
and shows a gradual shift of the peak in the co-moving 
space density to lower redshifts with declining luminosity. 
This behavior may be evidence for the growth of lower 
mass black holes emerging in an "anti-hierarchical" or 
"cosmic downsizi ng" fashion while accreting near their 
Eddington limit ([Marconi et al.ll2004l: fMerloni 20041. or 
the embers of a fully matured SMBH population with 
sub-Eddington accreti on rates. The former scenario has 
been substantiated bv IBarger et al.l ([20051 ) based on the 
optical lu minosities of the galaxie s hosting X-ray selected 
AGN, and lHeckman et all ([2004') using type 2 AGN from 
the SDSS with the [OHI] emission line luminosity as a 
proxy for the mass accretion rate and an estimate of the 
black hole mass from the M — a relation. 

Even though there has been much progress, there are 
remaining uncertainties in the current measure of the 
XLF. (1) A significant number of X-ray sources in the 
recent sur veys with Chandra an d XMM-Newton are not 
identified. iMainieri et al.l ([2005f ) find that the peak of the 
redshift distribution shifts to higher values (z ~ 1.2 — 1.5) 
when incorporating photometric red shifts for objects too 
faint for optical spectroscopy. (2) .Barger et al.l ([2005D 
demonstrate that the XLF can be fit equally well by a 
PLE model at z < 1.2. These models only begin to sub- 
stantially differ for low luminosity AGN at z > 1.5 where 
AGN statistics are quite low with most being provided 
by the deep CDF-N and CDF-S observations. New mod- 
erate depth surveys suc h as the Exteiided C handra Deep 
Field-South (E-CDF-S: [Lehmer et al.|[2005l ) and the Ex- 
tended Groth Strip fEGS: iNandra et all 12005) wiU pro- 
vide additional AGN at these luminosities and redshifts 
but await optical foUowup. (3) How does the AGN pop- 
ulation behave at redshifts above the peak (z ^ 2.5) 
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of the optically- selected QSO population? We have pre- 
sented evidence (|Silvemianll2004l : [Silverman et al.ll2005lJ ) 
for a similar evolution of luminous X-ray selected QSOs 
to th ose found in the optical surveys (e.g. .[Ri chards et al.l 
I2OO6I ) with a decline in the co-moving space density at 
z > 3 but these AGN are mainly type 1. We don't ex- 
pect the inclusion of luminous absorbed QSOs to dras- 
tically alter our measure of the XLF since t hey may be 
atmost as numerous a,s the type 1 QS Os (jGilli et al.1 
I2007D . Recent radio (IWall et al.l \2005) and near in- 
frared ([Brown et al.ll2006l ) survevsare further supporting 
a strong decline in the co-moving density at high redshift. 

In the present study, we measure the XLF of AGN in 
the hard X-ray (intrinsic 2-8 keV) band with an empha- 
sis on reducing uncertainties at high redshift (3 < z < 5). 
This paper is an extension to the preliminary results on 
the CO- moving s pace density of AGN as reported by th e 
ChaMP survey ([Silverman[[2004[ : [Silverman et al.[[2005bD . 
As previously described, these early epochs represent the 
emergence of the luminous QSOs and the rapid growth 
phase of young SMBHs. To date, the limited num- 
bers of X-ray selec ted AGN at z > 3 have constrained 
curre i it measures ([La Franca et al.l [20051 [Barger et al.l 
[2005t [Ueda et al.[ [2003[) to lower redshifts. Motivated 
bv [Barger et al.l ( 2005[ ). we use the observed soft X-ray 
band for AGN selection above z = 3 where we measure 
the rest-frame energies above 2 keV. Due to the rarity of 
luminous high redshift AGN, such an endeavor requires a 
compilation of surveys that covers a wide enough area at 
sufficient depths. As previously mentioned, a large dy- 
namic range from the deep, narrow pencil beam to the 
wide, shallow surveys is required to measure a luminos- 
ity function that spans low and high luminosities at a 
range of redshifts. Currently, there are a handful of deep 
surveys with Chandra (i.e., CDF-N, CDF-S) and XMM- 
Newton (Lockman Hole) that have published catalogs 
with a fair sample of low luminosity (42 < log Lx < 44) 
AGN out to ^ '--^ 5. To provide a significant sample of 
more luminous AGN {log Lx > 44), the Chandra Mul- 
tiwavelength Project (ChaMP) is carrying out a wider 
area survey of archived Chandra field s and the CLASXS 
([Yang et al.[[200l [Steffen et al.l[200l survey is imaging a 
contiguous area with nine Chandra pointings. The statis- 
tics of high redshift AGN are sure to improve with the 
anticipated results from the S WIRE/ Chandra (Wilk es et 
al., in pre paratio n), X Boot es ([Murray et al. [ 20051) . E - 
CDF-S (Lchmcr et al. 2005V EG S ([Nandra et al.[ [20051 ). 
XMM/COSMOS (Hasinger eraI1[2007f l. and the newly 
approved C/iandra/COSMOS surveys. 

We organize the paper as follows: Section |2l we de- 
scribe the various surveys used in this analysis including 
X-ray sensitivity, sky area coverage, incompleteness as a 
function of not only X-ray flux but optical magnitude, 
and AGN selection. Our method for measuring the lu- 
minosity function is presented in Section |3| and the re- 
sults, including best-fit analytic models, in Section |4| for 
all AGN types. In Section HI we address the resolved 
fraction of the CXRB and any underrepresented source 
populations. We directly compare our luminosity func- 
tion to that of optically-selected samples in Section |6l 
Based on our luminosity function, we derive in Section |7| 
the accretion rate distribution as a function of redshift 
and the cumulative mass density of SMBHs. Section [7Tl 
we compare the global mass accretion rate of SMBHs to 
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Fig. 1. — Observed flux in soft X-rays as a function of redshift 
for a theoretical AGN (F = 1.9) of luminosity (units of erg s~^) 
as labelled. We do not consider intrinsic absorption that may be 
present in many AGN. We highlight the redshift range 3 < ^ < 6 
over which an extension of the known XLF is needed. The hori- 
zontal lines mark the flux limits of various X-ray surveys capable 
of detecting AGN at these high redshifts. 

the star formation history of galaxies out to z 5. We 
end in Section [5] with some predictions of the numbers of 
high redshift (z > 3) AGN expected in new surveys that 
effectively enable us to extend the luminosity function 
to these high redshifts with accuracy. Throughout this 
work, we assume Ho — 70 km s~^ Mpc~^, = 0.7, and 
= 0.3. 

2. COMPILATION OF AGN FROM VARIOUS X-RAY 
SURVEYS 

The catalog of X-ray sources from the ChaMP pro- 
vides the foundation for our measure of the AGN lumi- 
nosity function. To this, we incorporate available cata- 
logs from additional X-ray surveys, described below and 
listed in Tabled! that effectively improve the coverage of 
the luminosity — redshift plane. In Figure (U we illus- 
trate the complementarity of surveys of various depths 
to probe the high and low luminosity ends of the popula- 
tion; the soft band is shown to highlight the feasibility of 
these surveys to detect AGN out to high redshifts {z ~ 5) 
by taking advantage of the higher sensitivity of the soft 
band compared to harder {E > 2 keV) energy bands. As 
is evident, the deep surveys (i.e., CDF-N, CDF-S) are key 
to improving the quality of our XLF at low luminosities 
(e.g., log Lx ~ 42 at 1 < 2; < 3) thus constraining a pre- 
viously unexplored part of the luminosity function before 
the Chandra and XMM-Newton era. The ChaMP and 
CLASXS surveys, of shallower depth but of wider area, 
effectively supply the more luminous {log Lx > 44) AGN 
especially at z > 1.5 that are underrepresented in the 
aforementioned deep fields due to a steeply declining lu- 
minosity function. These wider-area surveys are required 
since a sky coverage of over one square degree is needed 
to detect significant numbers of these rare, high redshift 
AGN. 

We generate both a hard (2-8 keV) and soft (0.5-2.0 
keV) X-ray se lected, flux-limit ed source catalog as mo- 
tivated by Ba rger et al.[ ([2005[ ) to determine their rest- 
frame 2-8 keV luminosity for sources with available red- 
shifts. We thereby optimize uniformity in regards to se- 
lection across the entire compilation of data sets. The 
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hard-band detections are used to construct a low-redshift 
sample at z < 3. This enables our selection to be less 
affected by absorbing columns up to Nh ~ lO^"^'^ cm~^. 
The soft band allows us to take advantage of Chandra^s 
high collecting area at low energies to detect faint, high- 
redshift AGN at z > 3. Above this redshift, intrinsic 
absorbing columns, if present, will have less effect on the 
observed flux since we are detecting X-rays at rest-frame 
energies above 2 keV. This selection technique also min- 
imizes uncertainties induced by large fc-corrections. 

Our sample has additional selection criteria due to the 
dependence of source identification on optical magnitude. 
This is clearly evident in the ChaMP by the fact that 
the identification of optical counterparts is currently re- 
stricted to i' < 23.5 (|breen et al.ll2004l : [Silverman et al.l 
l2005bD . This is more severe {i' < 22) when considering 
those optical counterparts that have reliable redshifts. 
Additional surveys included in this study are required to 
have optical imaging that is available that covers a wave- 
length range similar to the SDSS r' and i' bands used in 
the ChaMP. For hard-selected AGN < 3), the r' band 
is sufficient to measure their optical brightness. We use 
the i' band for the soft-selected sample {z > 3) since 
at these high redshifts extinction due to the intervening 
IGM becomes significant for bandpasses at lower wave- 
lengths. The addition of the Chandra Deep fields not 
only pushes our overall catalog to fainter flux limits but 
improves our characterization of X-ray bright, optically 
faint sources not yet identified in the ChaMP. 

Our selection method based on both the X-ray and 
optical flux is illustrated in Figure [H The fraction of X- 
ray sources with redshifts, either spectroscopic or pho- 
tometric, is represented as a grey scale image deter- 
mi ned using the a d aptive binning procedure described 
in [Silverman et al.l ()200 5tf) , and a function of both X- 
ray flux and optical magnitude. For the hard selected 
sample {top panel), we include X-ray sources from all 
surveys with /2_8 keV > 2.7 x 10~^^ erg s~^ cm~^ and 
r' < 22.0. To this, we add sources with fainter X-ray 
fluxes (/2-8 kcv < 2.7 X 10~^^ erg s~^ cm~^) and op- 
tical magnitudes (r' > 22.0) from all surveys with the 
exception of the ChaMP and AMSS. We set X-ray flux 
(/2-8 kov > 6.3 X 10~^^ erg s~^ cm~^) and optical mag- 
nitude limits (r' < 24.0) that effectively exclude the 
zone enclosed by the solid contour, where < 50% of the 
sources have redshifts. The same sort of two zone selec- 
tion [bottom panel) in the fx / opt plane occurs for the 
soft-selected sample with the inclusion of all sources hav- 
ing /0.5-2 koV > 1-0 X 10"i5 erg s"^ cm"^ and i' < 22.0; 
sources with fainter X-ray fluxes and optical magnitudes 
do not include ChaMP detections since currently few 
have identifications. For the soft-selected sample, we use 
the i' magnitudes to mitigate any flux loss due to extinc- 
tion blueward of Lya in high-redshift AGN. Again, we 
set X-ray flux (/o. 5-2.0 kcV > 1-0 x 10"^^ erg s"^ cm~^) 
and optical magnitude limits (i' < 24.0) where >50% of 
the sources in the fx /opt plane have redshifts. The flux 
distribution of X-ray sources from our above selection is 
shown in both the hard (Figure ^ and soft (Figure d]) 
bands. We also show the distribution of those with op- 
tical counterparts and identiflcation (i.e., redshifts) with 
spectroscopic or photometric techniques. Only the CDF- 
N and CDF-S surveys include photometric redshifts due 
to the rich data sets in these fields. As is evident, we 
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Fig. 2. — Fraction of identified (i.e., redshift is known) sources 
as a function of X-ray and optical flux in the hard/r' (top) and 
soft/i' (bottom) energy bands. The grey scale image ranges from 
zero (darkest) to 100% (lightest) of sources having redshifts mainly 
from optical spectroscopy. The AGN are shown with a black dot 
in the hard map for those at 2: < 3 and soft map for those at 2: > 3. 
In the soft band map, those AGN at z > A are also marked by 
a square. The X-ray and optical flux limits for our analysis are 
shown by the dashed and dotted lines. Since the ChaMP survey 
has a reasonable degree of completeness only at brighter flux levels, 
we only include ChaMP AGN that fall in the lower right region as 
defined by the dashed lines. The faintest limits (dotted line) are 
determined by the fraction of identified sources (> 50%). 

have redshifts for over ^ 50% of the X-ray sources at all 
fluxes. 

High X-ray luminosity becomes our single discriminant 
for AGN activity since many of the traditional optical 
AGN signatures are not present in obscured sources or 
may be outside spectroscopic coverage. The rest-frame 
2.0-8.0 keV luminosity (uncorrected for intrinsic absorp- 
tion) is calculated by assuming a fc-correction based on 
a powerlaw X-ray spectrum with photon index F = 1.9 
and required to exceed 10^^ erg thereby excluding 
any sources with X-ray emission primarily from a stellar 
or hot ISM component. This fixed lower limit on the 
luminosity {L2-S kcV > 10^^ erg s~^) is uniformly ap- 
plied to all X-ray sources from surveys considered and 
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Fig. 3. — top X-ray flux distribution for the sources selected 
in the 2-8 keV band: all X-ray sources (solid line), those with 
optical counterparts (r'; dashed line), optically-faint counterparts 
only (r' > 24; dash-dotted), and identifications based on a spec- 
troscopic or photometric redshift (dotted line), bottom Optical 
magnitude distribution. In both plots, the vertical lines show our 
chosen flux/magnitude limits with the brighter one marking the 
self-imposed limit for the ChaMP survey. 

yields a sample of 682 AGN with the contribution from 
each survey described below and listed in Table [TJ In 
Figure [HI we show the luminosity-redshift distribution 
of the full sample. We now have a significant sample of 
31 AGN at z > 3 to evaluate the luminosity function at 
high redshift. For the subsequent analysis, we make no 
attempt to subdivide the sample into the usual optical 
or X-ray based categories of AGN (e.g., type 1/type 2, 
obscured/unobscured). We do show for illustrative pur- 
poses in Figure El those AGN that have been classified 
through spectroscopy or photometry as a type 1 AGN 
(filled symbol). It is worth highlighting the fact that 
almost all of the AGN at z > 3 may be optically unob- 
scured as characterized by the presence of broad emission 
lines (FWHM > 1000 km s~^) in their optical spectra. 
Further spectroscopic foUowup with deep observations on 
8m class telescopes are required to adequately assess the 
contribution of obscured AGN at these redshifts. 

2.1. Chandra Multiwavelength Project (ChaMP) 

The ChaMP* , a survey of serendiptious, extragalactic 
X-ray sources detected in fields found in the Chandra 

* http://hea-www.harvard.edu/CHAMP/ 
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Fig. 4. — Same as Figure [3] though for the soft 0.5-2.0 keV and 
i' selected sources. 

archive, covers a large enough area at flux depths re- 
quired to detect significant numbers of AGN at high red- 
shift. We refer the reader to the following ChaMP papers 



for full de tails c oncerning the X-ray analvsi s (Kim et al 



2007a, b, 2004allH) and optic al foUowup ([Green et al 
,2004: .Silverman et all l2005aD_ programs. Br iefly, the 



ChaMP point source catalog (|Kim et al.l[2007ah contains 
7,365 unique sources found in 149 ACIS fields with expo- 
sure times ranging from 0.9 to 124 ksec that corresponds 
to a flux limit of /o. 5-8.0 keV = 9 x 10~^^ erg cm~^ s~^. 
The area coverage reaches ~ 10 deg^ at the brightest 
fluxes. The backbone of the optical foUowup program is 
a deep imaging campaign with the CTIO Blanco 4m and 
KPNO Mayall 4m through NOAO providing images in 
the SDSS filters g', r' , and i' with limiting magnitudes 
(> 24) depending on the depth of the X-ray exposure. 
We have undertaken an optical spectroscopic foUowup 
program utilizing the facilities available through NOAO 
(i.e., WIYN/3.5m, Blanco 4m, Gemini-N) and SAO (i.e., 
Magellan, MMT, FAST/1. 2m). Recently, we have be- 
gun a deeper spectroscopic campaign on Magellan with 
longer exposure times (~ 4 hours) and Gemini-N obser- 
vations to identify the optically faint counterparts that 
should include the type 2 QSOs that are expected to 
be found in ChaMP in g reater numbers (see Figure 13 
of [Silverman et al.ll2005a[ ) and the most distant (z > 4) 
QSOs. In total, we have accumulated ~ 1000 secure 
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Fig. 5. — Luminosity-redshift distribution of AGN from the X-ray surveys as labelled. Of the 682 AGN, 414 have broad emission lines in 
their optical spectra shown by the solid pentagon. Open squares correspond to AGN characterized by narrow optical emission or absorption 
lines. The white area highlights the regions for which our ID fraction is sufficient to measure the density of AGN. The dashed lines denote 
the Lx — z bins used for the 1/Va method. 



spectroscopic redshifts of X-ray sources to date. In the 
near future, the ChaMP anticipates the release of the op- 
tical imaging and photometry for 66 fields (W. Barkhousc 
et al., in preparation) and spectroscopy products (Green 
et al., in preparation) to the community. The ChaMP is 
also now extending to a total of 392 fields, using Chandra 
Cycle 2-6 observations overlapping the SDSS (Green et 
al. 2008). We have also begun to acquire near-infrared 
observations in the J and Ks bands for select fields using 
ISPI on the CTIO BLANCO 4m to measure photomet- 
ric redshifts and near-infrared properties for sources not 
accessible to optical spectroscopy. 

Based on a comprehensive catalog of 1,940 X-ray 
sources detected in 25 ChaMP fields (Tabled]), we have 
selected a subsample of hard (2.0-8.0 keV) and soft 
(0.5-2.0 keV) sources subject to X-ray and optical se- 
lection. Fields^ were chosen to have a limiting flux ca- 
pable of detecting high redshift AGN, quality optical 
imaging in the r' band for the z < 3 sample and i' for 
the z > 3 sample, and a substantial amount of opti- 
cal spectroscopic followup. Sources with greater than 10 
counts in their respective selection band are included. 
We limit our sources selected in the hard band to 392 
with /2_8 kcv > 2.7 X 10"^^ erg cm^^ s~\ r' < 22.0 and 
the soft sample to 609 with /o. 5-2.0 kcV > 1.0 x 10~^^ erg 
cm~^ s""'^, i' < 22.0. These flux limits reflect the X-ray 
and optical parameter space for which we have identified 
the majority (~ 74%) of the X-ray sources through our 

^ Two of the twenty-five fields lack i' imaging and are therefore 
not included in the soft-selected sample. 



optical spectroscopic observations. The brightest opti- 
cal counterparts were targeted for spectroscopy. This is 
clearly evident in Figure[6]that shows that most (87%) of 
the optically- bright (< 22 mag) counterparts to both soft 
and hard sources have measured redshifts. In Table [31 we 
give the number of sources for which we have obtained 
optical spectra and those with reliable redshifts (quality 
flag greater than 1). A confidence level is assigned to 
each redshift: l=uncertain, 2=probably correct, 3=se- 
cure. Spectra with an assigned level of one tend to in- 
clude BLAGN with a single low signal-to- noise (S/N) 
broad emission line, whereas, those with higher S/N can 
be assigned a two or three depending on the number of 
visible spectral features. Details of our spectroscopic pro- 
gram including redshift measurements and quality assess- 
ment will be presented in a future ChaMP paper (Green 
et al. 2008). Spectroscopic redshifts from the ChaMP 
are sup plemented with an additional ~ 20 from the 
SEXSI (|Eckart et al.]l200l and CYDER (|Treister et all 
[2001 surveys that have carried out observations of X-ray 
sources in a subset of the ChaMP fields. We have further 
removed 5 objects identified as clusters based on th eir ex- 
tended X-ray emission fsee iBarkhouse et"aLll2006[) . 



The sample of 286 AGN with L2-8 keV > 10*^ 



erg 



s~^ from the ChaMP improves the statistics in select 
regions of the Lx — z plane and increases the number of 
rare types as follows: (1) type 1 or broad emission line 
AGN (BLAGN thereafter) at z > 1.5 detected in the 
2-8 keV band, (2) high redshift AGN, twenty-one AGN 
have been identified at z > 3 (4 at z > 4) in the entire 
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The on-axis flux limits of this 
.5 X 10^^^ erg cm^^ s^^ and 
.f o.5-2 keV = 3 6 X 10 "^^ erg cm-^ s"^ 

iBarger et al.l ()2003D have identified counterparts for 
all of the 503 X-ray sources with deep optical imaging 
( B .V .Rn .In .z' ■. \Ca.TDak et al.ll2004f). We use photom etric 
transformations ( |Fukugita et al.lll996l : ISmith et al.1l2002f) 
to determine r' and i' magnitudes in the SDSS photo- 
metric system for optical counterparts with V, Rc and 
Ic magnitudes. Using the 486 counterparts with mag- 
nitude in these three bands, we measure the following 
mean values: < r' - Rc >= 0.27, < i' - Ic >= 0.45, and 
< Rc — Ic >= 0.80. We use these to determine r' and 
i' when a either V, Rc or Ic is undetermined. Only one 
X-ray source has no detections in all three bands {V, Rc 
and Ic)- Any uncertainty associated with these trans- 
formations does not significantly impact our measure of 
the luminosity function since we mainly consider these 
magnitudes to measure the fraction of identified X-ray 
sources using wide magnitude bins. Redshifts are avail- 
able for 319 (74%) of the soft-band and 223 (68%) of 
the hard-band detections with the majority (~ 80%) of 
these obtained from optical spectroscopy. The remain- 
ing identifications are derived from photometric tech- 
niques to obt ain redshifts for the optically-faint sources. 
iBarger et al.l ([2003) show that these are primarily at 
z > 1 and have optical/near IR emission representative 
of early-type galaxies. Above our self-imposed X-ray flux 
and optical magnitude limits (Table [T]), the CDF-N sam- 
ple provides 112 AGN with L(2.o-8.o kcV) > 10"*^ erg 
of which 8 lie at z > 3. 



Fig. 6. — Optical magnitude distribution (solid line) of 793 X- 
ray selected sources from the ChaMP in the hard sample (r'; top 
panel) and 1,125 in the soft sample («'; bottom panel). In both 
panels, we show those sources that have available optical spectra 
(dashed line) and those with reliable redshifts (dotted line). See 
Table [3] for a full listing of the ChaMP number statistics. 

ChaMP program to date, roughly half the published X- 
ray selected high redshift population. Optical spectra are 
shown in Figure [7] for 8 of the 13 z > 3 AGN included in 
this study and no t yet published) , and (3) ^ 20 galaxies 
with log Lx ~ 43 (|Kim et al.ll2006[ ) and no strong optical 
emission lines (Equivalent width; W\ < 5 A). 

2.2. Chandra Deep Field North (CDF-N) 

The CDF-N, with an on-axis exposure time of « 
2 Ms, is the deepe st X-ray image ever obtained 
(| Alexander et al.ll2003l) . Over a narrow field (0.12 deg^), 
a sample of 503 X-ray sources are detected down to a 
limiting flux of « 2.5 x 10~^^ erg cm~^ in the soft 
band (0.5-2.0 keV) and « 1.4 x lO'^^ erg cm'^ s^^ in 
the hard band (2.0-8.0 keV). This is 8x fainter than 
the limiting flux of the deepest fields in the ChaMP sur- 
vey. With follow-up observations in the radio through 
ultraviolet, this combined dataset is one of the richest 
available for studies of the X-ray emitting extragalactic 
source population. 

We select sources from the X-ray catalog 
(|Alexander et al.l l2003f ) that have greater than 10 
counts in either band to exclude those with less certain 
count rates. This restricts the X-ray catalog to 330 
sources in the hard band and 433 in the soft band. 
All fluxes are corrected for Galactic absorption using 



2.3. Chandra Deep Field South (CDF-S) 

We incorporate X-ray sources detected in the 1 Msec 
Chan dra observation of the CDF-S (jCiacconi et al.l 
l2002t l with greater than 10 counts in either the soft or 
hard X-ray band. Sixteen sources were removed from the 
main catalog that had more than one possible optical 
counterpart and four extended sources associated with 
an optical group/cluster. This results in a sample of 222 
hard and 256 soft-selecte d X-ra y sources. The fluxes re- 
ported in lGiacconi et al.l ()2002f ) are corrected for a neg- 
ligible amount of Galactic absorption. The flux limits 



erg cm 



of this sample are /2-8 keV = 4.6 x 10~^^ 
and /0.5-2 keV = 5.6 x 10"^'' erg cm~^ s~^. 

We utilize the wealth of optical imaging taken of the 
CDF-S to determine effective SDSS magnitudes (i.e., r', 
i') , ob j ect t ypes and redshifts. The R magnitude, given 
in Giacconi et al. (2002), and the relation (< r' — Rc >= 
0.27) found using the CDF-N sources, are used to deter- 
mine r' . We use the positions of the optical co u nterp arts 
of the CDF-S sources given in IZheng et all (|2004D to 
search for HST/ACS i' counterparts in the GOODS data 
(Giavalisco et al. 2004). For sources outsi de the GOODS 
cover age, we use the ESO imaging survey (jArnouts et al.l 
I2001h to determine the / magnitude . In to tal, 301 (88%) 
X-ray sources listed in IZheng et al.l ()2004l ) have an opti- 
cal counterpart in the iACS or Ieso band. We convert 
the Ieso to i using < i' — Ieso >= 0.49, as measured 
from 148 sources with detections in both the GOODS and 
ESO imaging surveys. The ESO imaging survey provides 

10 PIMMS is mostly written and maintained by Koii Mukai 
i jhttp: / /heasarc. nasa.gov/docs /software / tools/ p imms.html} 
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Fig. 7. — Optical spectra of AGN at 2 > 3 identified in the ChaM P. These comprise 8 of the 13 high redshift AGN from ChaMP included 
in our analysis. The remaining four have been published elsewhere ISilverman et aLll2002l . [20031) . 



/ magnitudes for 100 of the X-ray sources not covered by 
the GOODS data. 

Rcdshifts are available for 207 of the 256 soft-band 
sources with many (115) from spectroscopic foUowup 
(jSzokoly et al. 2004i). Multi-band p hotometry confirm s 
the redshifts of many of these sources (jZheng et al.ll2004f ). 
including a f ew with uncertain spectroscopic redshifts 
(0 < Q < 2; ISzokolv et al.ll2004D . and identifies an ad- 
ditional 49 soft-selected X-ray sources. We incorporate 



the classifications using photome tric technique s in ca ses 
where the quality as defined in I Zheng et al.l ()2004f ) is 
0.5 < Q < 1.0. This range of quality defines those pho- 
tometric redshift estimates that have similar results from 
at least two of the three methods. We also include photo- 
metric redshifts for thos e sources from the COMBO- 17 
survey (|Wolf et al.l l2004f) with i? < 24 since these red- 
shifts are fairly accurate at these brighter magnitudes. 
The CDF-S contributes 97 AGN (2 at z > 3) with X- 
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ray fluxes and optical magnitudes satisfying our chosen 
limits (Table [J). 

2.4. Chandra Large Area Synoptic X-ray Survey 
(CLASXS) 

The CLASXS survey is a contiguous mosaic of 9 
Chandra ACIS-I pointings in the low Galactic column 
region o f the Lockman H ole-Northwest. The X-ray 
catalog (jYang et al.l l2004f ) contains 519 X-ray sources 
detected in the 0.4-2.0 keV band with a flux above 
5 X 10~^^ erg cm~^ s~^. To merge the sample with the 
ChaMP, we use the Ga l actic c orrected 0.5-2.0 keV fluxes 
listed in iSteffen et all (|2004f ) for sources with greater 
than or equal to 5 counts in the respective energy se- 
lection band. We note that the conversion from counts 
to flux takes into account the hardness ratio of each in- 
dividual source, similar to the CDF-N data. This can 
cause slight non-uniformity in our overall compilation. 
We assume that t he area coverage shown in Figure 9 of 
lYang et al.l (|2004j ) in the 0.4-2.0 keV does not change 
significantly when using a slightly wider energ y band. 
With deep optical imaging {B,V,Rc,Ic:z'), Steffe n et al.l 
(j2004l) have identified optical counterparts for 99% of 
the X-ray sources. As done for the CDF-N (Section [^21) 
sample, we use the photometric transformations to de- 
termine r' and i' magnitudes in the SDSS photomet- 
ric system for optical counterparts with V, Rc and 
Ic magnitudes. Using the 451 counterparts, we mea- 
sure: < r' - Rc >= 0.23, < i' - Ic >= 0.43, and 
< Rc — Ic >= 0.37. Only five X-ray sources have no 
detections in all three V, Rc and Ic ban ds. Optical spec- 
troscopic classifications are available (St effen et al.ll200^ 
for 52%, including a significant number of optically-faint 
(i? > 24) counterparts. The spectroscopic classification 
scheme is nearly identical to the ChaMP. We group the 
emission line objects (i.e., star forming and Seyfert 2 
galaxies) and absorption line galaxies under a general 
"galaxy" category. Above the flux limits (Table [T|) im- 
plemented for this study, there are 106 sources attributed 
to AGN activity based solely on their X-ray luminosity. 

2.5. XMM-Newton Lockman Hole 

We include X-ray sources found in t he deep XMM- 
Newton observation of the Lockman Hole (|Hasinger et al.l 
[2001.) to boost our sample at high redshift. We only 
use the soft (0.5-2.0 keV) selected sources to add to 
our z > 3 sample since the optical foUowup has pri- 
marily targeted sof t-selected sources from the ROSAT 
Ultra Deep Survey (jLehmann et al.|[2001h . The most re- 
cent catalog (Brunner et al. in preparation) contains 
340 sources detected in the soft band with a limiting 
flux of 2 X 10~^^ erg cm~^ s~^. We limit the X-ray de- 
tections to those within 12' of the aimpoint since this 
region contains the majority of the sources with spec- 
troscopic redshifts. This sample provides a fair number 
of AGN at high redsh ift since four have been reported 
(jLehmann et al.l[2001l ) at z > 3 and one additional AGN 
has been identified at z = 3.244 (Szokoly et al. in prepa- 
ration). Optical magnitudes have been converted to 
SDSS magnitudes as described above (Section 12. 2p . We 
further limit the catalog by selecting only those sources 
with X-ray flux above fx > 1.0 x 10^^'^ erg cm~^ s~^ 
since the fraction of identifled sources at fainter fluxes 



is low (< 30%). Above these flux limits, there is a rea- 
sonably high degree of completeness with 59% (69 out 
of 117) of the sources identified by optical spectroscopy. 
This field adds 5 AGN at z > 3 that also satisfies our 
optical magnitude selection {i' < 24). 

2.6. ASCA Medium Sensitivity Survey (AMSS) 

We include AGN from the AMSS (lAkivama et all 
l2003f ) to tie down the bright end of the luminosity 
function at low redshifts. This survey presents 87 
sources detected in the 2-10 keV band with fluxes 
down to 3 X 10"^'^ erg cm~^ s^^. We convert flux to 
that in the 2-8 keV band using PIMMS. We fur- 
ther restrict the sample to 76 AGN with fluxes above 
2.5 X 10""'^^ erg s~^ cm~^; t his is essentially the same 
limit chosen by lUeda et all ()2003f ) to mitigate the influ- 
ence of systematic effects near the flux limit of the survey. 
Optical photometry and source classifications are avail- 
able to easily assimilate this sample of 76 AGN into our 
composite catalog. 

2.7. Sky coverage 

Each survey included in this study has published a 
measure of the sky area as a function of flux in both the 
soft and hard energy band s. To this, we ad d the con- 
tribution from t he ChaMP (iKim et al.ll2007allbl V As de- 
scribed briefly in [Silverman et alT ( 2005b[ ) and reiterated 
here, a series of simulations are performed to character- 
ize the sensitivity, completeness, and sky area coverage 
as a function of X-ray flux. The simulations consists of 
three parts, 1) generating artiflcial X-ray sources with 
MARX (MARX Technical Manual") and adding them 
to real X-ray images, 2) detecting these artificial sources 
by using the CIAO3.0^'^ wavdetect and extracting source 
properties with XPIPE identically as performed for ac- 
tual sources, and 3) estimating the sky area coverage as 
a function of fiux by comparing the input and detected 
source properties. The simulations are restricted to spe- 
cific CCDs for ACIS-I (10, II, 12, and 13) and ACIS-S 
(12, 13, S2, and S3). These CCDs are closest to the aim- 
point for each observation. ACIS S4 is excluded because 
of its high background and noise streaks. Sources far 
off-axis (9 > 12') are excluded since the flux sensitiv- 
ity is low and the PSF is degraded. These simulations 
allow determination of corrections for the source detec- 
tion incompleteness at faint flux l evels quantified in the 
first ChaMP X-ray analysis paper (|Kim et al.ll2004bD . In 
Figure [H we show the sky coverage from all surveys in 
the soft (0.5-2.0 keV) and hard (2.0-8.0 keV) bands. As 
shown, the ChaMP adds ^ 2 square degrees at interme- 
diate flux levels between the narrow deep fields and the 
shallower, wide area ASCA surveys. 

3. MEASURING THE HARD (2-8 KEV) XLF 

We measure the differential XLF (d<I>/dlogL) expressed 
in Equation [T] where N is the number of AGN per unit 
co-moving volume (V) and log L^ as a function of X-ray 
luminosity L^ and redshift z. 



d^TV 



d^{L^,z) ^ 

dlogLx dV AlogL- 



(1) 



1^ http://space.mit.edu/CXC/MA RX 
|http: //cxc .harvard. edu/clao] 
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Fig. 8. — Sky area coverage as a function of X-ray flux in the 
hard (2.0-8.0 keV) band top and soft band (0.5-2.0 keV) bottom. 
The thick line is the sum of aU surveys. 



This function is assumed to be continuous over the range 
of luminosity and redshift spanned by our sample. The 
differential luminosity (dlogL) is expressed as a loga- 
rithm (base 10) due to the 4 orders of magnitude covered 
by our sample. 

Many methods to determine the luminosity function 
of extragalactic pop ulations have been implemented (see 
iBinggeli et al.llT988l . for a review) . The method of choice 
can depend on the complexity of the selection function 
and/or completeness of the sample. For example, the 
traditional l/Va method applied to a sample of quasars 
((Schmidt 1968) easily allows one to incorporate bivari- 
ate selection, mainly driven by incomplete optical spec- 
troscopic identification of objects selected in another 
wavelength regime (i.e., radio, X- ray) . Alternative tech- 
niques usi ng binned data (e.g., iPage fc Carreral [19991 : 
iMivaji e"t~a l. 2001) have been applied to X-ray selected 
AGN surveys, though they require fairly complete (i.e., 
optical identificatio n and redshifts) samp les. Maximum 
likelihood methods ([Marshall et al.|[T983( ). though model 
dependent, are most applicable when incorporating com- 
plex selection functions. They also mitigate effects that 
result from finite bin widths and irregular sampling of 
each bin, typically wors t near the flux limit (see Figure 
1 of iMivaii et al.l[200lh . 

We measure the luminosity function using two methods 
to provide further assurance of our results: (1) the 1/Va 
method and (2) a model-dependent maximum likelihood 
technique. The latter allows us to compar e model fits 
and best-fit parameters to those published bv lUeda et al.l 



(|2003f) . lBarger et al.l (|2005f) . and lLa Franca et all (|2005f) . 
In this work, we do not consider the absorption of X- 
rays intrinsic to our sample thus making no attempt to 
determine an Nh distribution as done i n recent studies 
(|Ueda et al.ll2003l : iLa Franca et al.ll2005f ) since it is diffi- 
cult to accurately measure absorption in AGN at these 
high redshifts {z > 3) due to the very limited count 
statistics. This choice does not severely affect our re- 
sults since we are using the 2-8 keV energy band and 
at these high redshifts we mainly probe the bright end 
of the luminosity function {log Lx > 44) that has been 
found to have a low fraction of ab sorbed sources (~ 34%; 
see Figure 7 of|UedaetaD[2Q0l). 

3.1. Binned 1/Va method 

We estimate the XLF in fixed lumi nosity and red- 
shift bins using the 1/Va method (jSchmidtl 119681 : 
lAvni fc Bahcalllll98Cll ) and our sample described in Sec- 
tion [51 For each L — z bin marked in Figure [5l the value 
of the XLF is a sum (Equation ^ of the contribution 
from each AGN falling within this specific bin. 



d<^>{L^,z) 
dlog Lx 



1 



N 



AlogLx j-[ 



(2) 



We exclude the low luminosity bins in each redshift in- 
terval that are not well sampled at our flux limit. Our 
sample at z > 3 is now sufficient to measure the lumi- 
nosity function in two redshift intervals. This enables us 
to further probe SMBH growth at these early epochs. 

We apply a correction factor (Ci) for each AGN to ac- 
count for the incompleteness in our optical spectroscopic 
identifications. This factor is the reciprocal of the frac- 
tion of identified sources (/id) at X-ray fiuxes and optical 
magnitudes comparable to each source, as shown by the 
grey scale image in Figure [2j The accessible volume Va 
is a f unction of bo th X-ray and optical limiting fluxes. 
From iHogd ([1999f ). we measure the co-moving volume 
for each AGN using the concordant cosmological model 
for nk = 0: 



C 



(»A,(l-Hz)3 + r!A)-V2 

(1 + Z)2 



(3) 



X V,a{Lx,i,z)dz 

The solid angle Oa(Lx,i, z) depends on the observed flux, 
that varies across the width (dz) of the redshift bin, of 
each AGN with specific X-ray luminosity {Lx,i) to deter- 
mine an accessible volume. This method does not con- 
sider evolution across the bin that may not be negligible 
when using large luminosity and redshift bins. Since we 
have a sample selected in two different energy bands, the 
sky coverage for AGN at 2 < 3 depends on their 2-8 keV 
X-ray flux as shown in Figure [5| top and their 0.5-2.0 
keV flux for objects at higher redshifts (z > 3; Figure [8] 
botto m). Further details o n our measure of /id is given 
in Sil verman et al.l ([2005bl ). We estimate Icr errors based 
on a Poisson distribution due to the small number of ob- 
jects per redshift bin. 

3.2. Analytic model fitting 

It is highly desirable to express the luminosity function 
as a smooth, well behaved analytic function in order to 
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maximize its versatility for a wide range of uses as later 
illustrated. A well defined model should be described by 
parameters that provide insight into the physical charac- 
teristics of the population (e.g., evolution rates). We fit 
our data with three different functional forms, as detailed 
below, over the full redshift (0.2 < z < 5.5) and lumi- 
nosity (42 < log Lx < 46) range spanned by our sample. 
We choose models that are conventionally used in the 
literature in order to directly compare best fit parame- 
ters and to extend their applicability to higher redshifts 
than have yet been accurately done in the hard X-ray 
band. As is usually the case with a data set that covers 
different parameter space, slight modifications of these 
well utilized models are required. 

To determine the values of the best-fit parameters for 
all mo dels, we implement a maximum likelihood tech- 
nique (|Marshall et al.lll983D . We minimize the fol lowing 
expre ssion using the MINUIT software package (jjamesl 
[l99l available from the CERN Pro gram Library. 



N 



i=l 



dlog ix 



(4) 



+2 1 I ^ Q{.fx,r'/i')— dz dlog Lx 



1 Jh dlog Lx 



dz 



Here, &{fx) is a correction for incomplete redshift in- 
formation as a function of X-ray flux as shown in the 
top panel of Figure [3] and |4l For this exercise, we have 
chosen to neglect the dependence of incompleteness as 
a function of optical magnitude. After each call to MI- 
NUIT the model is evaluated and compared to our total 
AGN sample size to determine the normalization {Ao). 
We estimate the 68% confidence region for each parame- 
ter while allowing the other parameters to float freely by 
varying the parameter of interest until A5(= Ax^) = 1.0 
(jLampton et al.lfr976( l. We note that this procedure may 
not accurately represent the true error interval since we 
are simultaneously fitting multiple parameters. For the 
normalization (Ao), we use the Poisson error based on 

the sample size instead of the value of AS. 

As shown by many studies (e.g., iMarshalll Il987t 
iBovle et"aI]|1988D , the shape of the XLF is best described 
as a double powerlaw modified by a factor(s) for evolu- 
tion. We first implement a 'pure' luminosity evolution 
model (PLE; Equation s [H - [7t t hat is almost id entical 
to that used by Richards et all (|2006l ) and Wolf et all 
(j2003). We drop the third order term for L* (Equa- 
tion [5]) because any incorporation of higher order evo- 
lution terms would require a significantly larger sample 
of X-ray selected AGN not yet available. In total, the 
model has 7 free parameters (Aq, 7I, 72, Lq, el, e2, Zc). 
Best-fit results for this model (PLE Model A) are shown 
in Table m from our maximum likelihood (ML) routine. 



d$(£x,^) ^ 

dLogLx iLx/L4z))'>^ + (Lx/L*(^))^' ' 



with 



logL* ^ logLo + elS_ + e2f , 



(5) 



(6) 



£, = log 



l + z 



(7) 



We are motivated to use a model that has more flex- 
ibility to fit our data, especially at low Lx- A modi- 
fied 'PLE' mo d el (E quations [5H7]) recently presented in 
iHopkins et al.l (|2007f ) is quite appropriate. This model 
adds two additional free parameters, a redshift depen- 
dence to the spectral indices of the double powerlaw. 
In our case, we choose to incorporate this feature only 
for the low luminosity slope (Equation [8|) that evidently 
changes with redshift. Best-fit results for this model 
(mod-PLE/Model D) are shown in Table [4] from our ML 
routine. 



72 = (72)0 



l + z 
1 + zc 



(8) 



Third, we implement a LDDE (Equations [5] ~[TTJ Note 
that for this case = Lq) model that has been use- 
ful in quantitatively describing the shape and evolution 
of the luminosity fun ction of X-ray selected AGN not 
only in the hard band (Ucda et al. 2003; La Franca et al.l 
2005f) but th e soft band (Hasingcr et al. 2005) as well. 
Mivaii et al.l (|2000D introduced a slight variant of this 
LDDE model to fit the AGN sample from a compilation 
of ROSAT surveys of varying depth and area coverage. 
Using this model, they were able to provide the first ev- 
idence that a luminosity-dependent evolution scheme is 
required for X-ray selected samples. 



d^Lx,z) ^ d^(Lx,0) 
dLogLx dLogix 



e(z, Lx) 



(9) 



where 



''^^'^"^"{i(i[(i + ^)/(i + ^c)]'^' \z>ll) ' 



(Z < Zc) 



and 



2c (ix) = 



_(z*{L^/L,r (ix<ia) 



(Lx > La) 



(11) 



In terms of model parameters, we aim to extend the ap- 
plicability of this model out to z 5 with a re-evaluation 
of the value of e2 that describes the evolution rate beyond 
the cutoff redshift Zc- From our prelimi nary measure of 
the cq -moving space density of AGN (jSilverman et al.l 
l2005bD . we e xpect this value to differ sharply from 
other studies (|Ueda et all 120031: iLa Franca et all l2005f) 
with a much stronger negative evolution rate that is 
similar to those found in optical QSO studies (e.g. , 
Warren et al. 1994:' Schmidt et al.lll995HWolf et al.l[2003l : 



Richards et al. 200(1 7" In addition, the values of a and 



z* may differ due to our larger AGN sample at high red- 
shift. To further si mplify this comple x model, we fix La 
to the value used bv lUeda et al.l (I2003D . In total there are 
8 free parameters(ylo,7l, 72, Lq, el, e2, 2*, a). Best-fit 
results for this model (LDDE Model B and C) are shown 
in Table m from our ML routine. 
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4. HARD XLF; RESULTS 
4.1. 1/Va method 

Our measure of the XLF, using methods detailed in 
the previous section, is presented in Figure [5] for seven 
redshift intervals with bins of finite width in luminos- 
ity. We show smooth analytic model curves (see fol- 
lowing section for details) including that evaluated at 
z = (dashed line) to aid in our visualization of the 
evolution with redshift and luminosity. First, we see 
that the shape of the XLF up to z 2 is well approx- 
imated by the familiar double powerlaw with a steeper 
bright end slope. The shape at higher redshifts is less 
constrained due to the limited statistics especially below 
the break or knee (i.e., characteristic luminosity dividing 
the faint and bright end slopes). To first order, there 
is an overall shift of the XLF to higher luminosities as 
a function of increasing redshift (i.e., pseudo-PLE) up 
to z ~ 3 as reported by many studies in various wave- 
bands over the past few decades. T here is a decline in 
either luminos i ty or space densitv (iBarger et all 120051 : 
iHasinger et"aLl l2005l : [Silverman et al.l l2005b[ ) at higher 
redshifts simila r to the beh avior s een in optical (e.g. , 
iFan et al.l [2OOII : IWolf et all ^ 2003: R ichards et all 12009) 
and radio-selected ([Wall et al.ii2005i ) samples. Our data 
preclude us from distinguishing these two modes of evo- 
lution at z > 3. Second, we see that there is evidence 
for a flattening of the faint-end slope with increasing red- 
shift most noticeable by comparing the data points with 
the z = model (dashed line) in the redshift interval 
1.0 < z < 1.5. The need to utilize more complex mod- 
els (i .e., LDDE), wh en dealing with X-ray selecte d sam- 
ples (IHasinger et al. 2005; Ucda et al. 2003; M ivaii et al.l 
|2000D . is primarily due to this behavior. 

In Figure [TU] we comp are our binned me asure of the 
XLF to the best-fit PLE (iBarger et al.ll2005D and LDDE 
(|Ueda et al.ll2003l : lLa Franca et al."20'05T models to check 
for consistency at z < 3. Barger et al. (2005) report 
that a PLE model best represents the data at z < 1.2, 
with characteristic parameters including e volution rates 
simil ar to optically-selected samples (e.g., ICroom et al.l 
|2004[) . As shown in the two low redshift panels {top; 
z < 1.5), it is apparent that the PLE model provides 
an adequate representation of the shape and evolution of 
the population. As previously mentioned, a flattening of 
the faint-end slope, most evident at 1.0 < z < 1.5 (top, 
middle panel) , is apparent and cannot be accounted with 
a PLE model. Our sample does not provide enough AGN 
with luminosities below the knee and z > 2 (top, right 
panel) to constrain the faint-end slope; this an unfor- 
tunate consequence of our optical magnitude limits and 
incomplete spectroscopic foUowup. In the lower panels, 
we compare our d ata with the more c omplex LDDE of 
lUeda et al.l (|2003D and iLa Franca et aL (2005). In gen- 
eral, the data agree well with both models that are nearly 
equivalent especially at high luminosities. It is worth not- 
ing that there is some slight uncertainty in the slope and 
normalization at the faint end {log Lx < 44) that may be 
due to X-ray sources that haven't been included due to 
our optical magnitude selection. We note that the lowest 
luminosity point in each of the panels is clearly affected 
by the fact that AGN do not entirely fill each Lx — z 
bin (see Figure O that falls on or near our flux limits. 
An accurate measure of the faint-end XLF is beyond the 



scope of this paper and requires more complete samples. 

As previously mentioned, we see a significant drop in 
the LF above z ~ 3 in either normalization or char- 
acteristic luminosity over two redshift intervals (Fig- 
ure [5]) . This behavior is similar to the decline seen in 
the soft and broad Cha ndra bands that we reported in 
ISilverman et al.1 (|2005bf) . With this new sample, we can 
further constrain the slope and normalization of the XLF 
at these high redshifts though still larger samples are re- 
quired especially at z > 4 where only 4 AGN are included 
in this analysis. Between 3 < z < 4, our sample of 26 
AGN does show a significant decline from the peak ac- 
tivity at 2 < z < 3. To check the integrity of our method 
that utilizes the observed soft X-ray band to measure the 
XLF above z = 3, we have measured the XLF in the red- 
shift interval 3 < z < 4 using the hard band as done at 
lower redshifts. These two measures, shown in FigurefTTl 
agree within ~ la for each data point. 

It is useful to illustrate how the ChaMP sample im- 
proves the measure of the hard XLF. In Figure [T^ we 
display the XLF as measured with the 1/Va method 
in the three highest redshift intervals. In the top row, 
the entire sample is used while the bottom row includes 
all AGN with the exception of those from ChaMP. The 
sample size is given for each data point. For re ference, 
the best-fit LDDE model of lU eda et all ()2003[ ) (sohd 
line) clearly demonstrates the difference in our results 
at z > 3. First, it is evident that a similar decline in the 
normalization of the XLF is seen with and without the 
ChaMP data at z > 3. The XLF in the bottom panels 
is in agreement with the behavior of the XLF reported 
by IBarger et al.l (|2005D using practically similar samples. 
The ChaMP does improve upon the accuracy of the high 
redshift XLF by boosting the overall numbers including 
those in specific regions of Lx — z plane having limited 
statistics: (1) the bright {log Lx > 45) end slope at 
2 < z < 3 is now further constrained, (2) measurement 
errors in the redshift interval 3 < z < 4 are effectively 
reduced, and (3) the numbers of z > 4 AGN are dou- 
bled with an additional data point at higher luminosities 
{log Lx ^ 44.8) that further constrains the slope. Re- 
sults ba sed solely on the Cha MP sample are presented 
in both ISilverman et al.l (|2005b ) and the PhD thesis of 
ISilvermanI (l200i ). 

As frequently represented, the XLF can also be plotted 
as a function of redshift for AGN within a fixed luminos- 
ity interval. We show in Figure [T3] the co-moving space 
density for three luminosity ranges. The lowest luminos- 
ity range is limited to z < 1 due to the lack of X-ray 
sensitivity at higher redshifts. For higher luminosities 
{log Lx > 43.5), our sample provides suitable numbers 
of AGN out to z ~ 5. We clearly see an increase in 
the number of AGN with redshift compared to the lo- 
cal universe followed by a decline beyond a peak whose 
redshift is dependent on luminosity. The most luminous 
AGN {Lx > 44.5) peak at a redshift around z ^ 2, sim- 
ilar to the optical surveys, while the lower luminosity 
AGN are most prevalent at z ~ 1.0. We have plotted 
the measurements using the observed hard band (open 
symbols) to show that the decline at z > 3 is not due to 
a significant population of high redshift sources not de- 
tected in the soft b and. These results are in agreement 
with similar stu dies (lUeda et al.ll2003l: IBarger et al.l l2005l: 
ILa Franca et al.l 120051 : IHasinger et al.ll2005D with" an im- 
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Fig. 9. — X-ray luminosity function. Tiic results from the binned 1/Va method are given by filled circles with Icr errors. Our best-fit analytic models using an unbinned maximum 
likelihood method are shown by the lines (solid=mod-PLE-Model D; dotted=LDDE-Model C). The dashed line is our LDDE model C evaluated at z = 0. 
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Fig. 10. — Comparison of our binned XLF (l/Va method; filled points) and analytic model fits at 2 < 3 to recently published models. 
Top ChaMP PLE model A (solid line) in three redshift inte rvals com pared to Barger ct al. (2005) (dot ted line) in the first two panels. 
Bottom Our best fit LDDE model B (solid line) compared to lUeda"eFal.i (|2003i ') (dashed-dotted line) and lLa Franca et al.l ((2005i) (Fit #4; 
dotted line). The dashed line is the Ueda model evaluated at 2: = in all panels. 



Evolution of X-ray-selected AGN 



15 



1 C'- 



1 0-5 



1 0- 



1 0"^ 



1 0- 



Soft selected: 6 
Hard selected: 2 




42 43 44 45 

Log [Lx (2.0-8.0 keV; erg s-')] 



46 



Fig. 11. — XLF in the redshift range 3 < 2; < 4 for AGN selected 
in the soft band (filled circles) and hard band (open circles). The 
numbers of AGN in each bin are reported at the top. The similarity 
of these XLFs and the improved statistics in the soft band illustrate 
our justification for using an energy dependent selection function . 
For reference, we show the analytic model from lUeda et all (pOO^) 
at 2 = (dashed line). 

provement in constraints at z > 3 for the hard XLF. 

We note that X-ray sources without redshifts con- 
tribute to the uncertainty in both our XLF and space 
density. While those which were never observed spectro- 
scopically might be expected to be similar in properties 
(e.g., in distributions of object type and redshift), the 
(significantly smaller) number of objects that we have 
observed spectroscopically without achieving a redshift 
may differ, which could bias our simple fractional incom- 
pleteness correction scheme. To illustrate their possible 
influence, we measure maximum values of the co-moving 
space density (Figure [T5)) by placin g all possible uniden - 
tified objects into each r edshift bin iBarger et al.l (|2005D : 
[Silverman et al.l ()2005bl ). We see that the highest red- 
shift bins (z > 3) are most susceptible to additional un- 
certainty if the unidentified objects have a redshift dis- 
tribution dissimilar to our AGN sample. 

4.2. Analytic model fitting 

Before we present a global fit over all redshifts and 
luminosities spanned by our sample, we verify that our 
method as described in Section [3T2l is robust. This is done 
by fitting our data with well known models over a similar 
range of redshift and then comparing the resulting best- 
fit parameters to published values. First, we ran our ML 
routine over a redshift interval of . 2 < z < 3.0 using a 
PLE model, similar to IBarger et al.l (|2005D and a LDDE 
model with fixed param eters (e2, Lg , and Zc) having the 
same value as given in lUeda et all (p003l ). Specifically 
for the LDDE model, we are interested in comparing the 
values of our free parameters descriptive of the shape 
(7I, 72), evolution (el) below the cutoff redshift {zc), 
and strength (a) of the dependence of Zc on luminosity. 

In Figure [TUl we plot the best-fit results from our ML 
routine (top^FLE Model A; 6ottoTO=LDDE Model B; 
see Table H] for actu al best-fit pa rameter s) and compare 
to published results (IBarger et al.,,2005. : La Franca et al.l 
120051 lUeda et al.ll2003l ) at z < 3 Three redshift inter- 



vals are shown in each case. Our P LE model (top panel) 
fit, in general, agrees with that of IBarger et all (|2005l ) 
at z < 1.5. Our fit at the faint end for z > 1 does not 
smoothly transverse the binned data points thus does not 
adequately represent the slope below the break. Over- 
all, there is a good agreement be tween our LDDE m odel 
(B; bottom panel) and that of lUeda et al] (|2003f) and 
iLa Franca et al.l (|2005( ) with some minor differences at 
the faint end that are worth noting. Below the knee of 
the luminosity function, our best-fit model falls below 
that of these two aforementioned measures. As previ- 
ously noted, this may be a result of our optical magnitude 
limits that remove a significant fraction of the hard and 
soft X-ray sources above our chosen X-ray limits. This 
may explain (1) our lower normalization (Ao), though 
sensitive to the values of the other parameters, as com- 
pared to both the Ueda and La Franca models, and (2) 
the differences between the faint-end slope at z > 1. The 
evolution factor (el) and power law index for the luminos- 
ity dependen t cutoff redshift (a ) are within the la errors 
reported bv lUeda et al.l ()2003f ). Given these small dis- 
crepancies and the limitations of the data, we conclude 
that all three models are equally valid at z < 3. 

We aim to extend the fit of the LDDE model to 
higher redshift (z > 3 ). This has been attempted by 
ILa Franca et al.l (j2005f ) though with a smaller sample 
at these redshifts that may also be biased by cosmic 
variance since more than half (5 of 9) of the high red- 
shift AGN are found in the limited area coverage of 
the CDF-S. Since the size our sample is also restric- 
tive at these redshifts (31), we need to reduce the num- 
ber of free parameters. To do so, we fix those that are 
highly constrained at lower redshift: 7I, 72, Lo, and La- 
The cutoff redshift (zc), evolution parameter above the 
cutoff redshift (e2) and a are free to vary. The most 
likely fit over the full redshift range 0.2 < z < 5.5 
is shown in Figure |2| (dotted line) with the values of 
each parameters given in Table |3| (Model C). First, 
there is agreement with the binned values from the 1/Va 
method with some exceptions elaborated further below. 
The main improvement, as compared to published hard 
XLFs, is that the fit requires a much stronger evolu- 
tion (e2 — — 3.271q 34) above t he cutoff redshift than 
previously found bv lUeda et al.l (j2003l e2 = —1.5) and 
iLa Franca et al.l ([20051, e2 = -1.15). Our value IS more 
constrained as evident by the la errors that are ~ 2 x 
smaller than those reported bv lLa Franca et alj (|2005f ). 

There is a slight discrepancy between the LDDE fit and 
the binned data from our 1/Va method worth highlight- 
ing. The strong evolution (e2) above the cutoff redshift 
from our fit causes the faint end slope {log Lx < 44.5) 
to flatten a bit more rapidly with increasing redshift 
than shown by the binned data. This is evident in the 
2 < z < 3 and 3 < z < 4 redshift bins in Figure H] with 
the lowest luminosity points falling above the LDDE fit 
(dotted line); although, the statistics are too low to def- 
initely make such claims. This may just illustrate the 
difficulties in accurately modelling the XLF with under- 
lying complexities in addition to the abrupt change from 
positive to negative evolution. It is also worth recogniz- 
ing that the faint-end slope of the luminosity function 
at these high redshifts in not well constrained since only 
the CDF-N and CDF-S are capable of detecting these 



16 



Silverman et al. 



Full sample 

6 25 30 5 




42 43 44 45 46 

Log Lx (2.0-8.0 keV) 



6 9 10 1 




42 43 44 45 

Log Lx (2.0-8.0 keV) 




42 45 44 45 

Log Lx (2.0-8.0 keV) 



Full sample without ChoMP 
5 17 13 1 




42 43 44 45 46 

Log Lx (2.0-8.0 keV) 




43 44 45 46 

Log Lx (2.0-8.0 keV) 




42 43 44 45 

Log Lx (2.0-8.0 keV) 



Fig. 12. — A comparison of the high redshift XLF between the full sample (top) and that ex c luding the ChaMP AGN (bottom). The 
data points are the same as in Figure [9] The solid line is the best- fit LDDE model of lUeda et al.l I I2003I) that exemplifies our differences at 
high redshift. The number of AGN per luminosity and redshift bin is shown above their respective data point. 
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Fig. 13. — Go-moving space density for three luminosity ranges using the 1/Va method. The open points show the measured values using 
the observed hard band. The lines (dotted and dashed) illustrate the maximum uncertainty by assuming that all possible unidentified 
objects fall in their respective bins. 
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AGN and their optical identification is quite challeng- 
ing. Some progress has been n iade using type I AG N se- 
lected in the soft X-ray band (|Hasinger et al.l [20051 ) and 
optically -selected QSOs fo und in Lyman break galaxy 
surveys (Hunt et al.l l2004[ ) that clearly show a flatter 
faint-end slope compared to lower redshifts. 

Finally, we use a modified PLE (mod-PLE) function 
as described above and in Table |4] to allow us greater 
flexibility. Accounting for the flattening of the faint- 
end slope (a) without a dependence on the evolution 
(e2) above the cutoff redshift may mitigate the problem 
mentioned above. The results shown in Figure [9] (solid 
line) are also consistent with our binned 1/Va measure- 
ments and show better agreement in the redshift interval 
2 < z < 3 than the LDDE fit. We find the dependence 
of the faint-end slope on redshift (a = — 1.04j^g |^2) a 
bit stronger than that measured bv I Hopkins et al.l ()2007L 
k^i = -0.623 ± 0.132). At higher redshifts (z > 3), this 
model may overestimate the numbers of low luminosity 
AGN log Lx < 44. This is evident by comparing the 
observed and expected numbers of high redshift AGN in 
the deep Chandra Fields (see Section [8]). Interestingly, 
the difference between these two models below the break 
luminosity becomes larger at z > 4 with ~ 4x more 
AGN with log Lx ~ 43. This region of Lx — z param- 
eter space is primarily an extrapolation of the data. In 
Section [HI we investigate whether it is possible to accu- 
rately measure the faint-end slope at these high redshifts 
with current or future Chandra observations. 

5. MISSED SOURCES, TYPE 2 QSOS AND THE COSMIC 
X-RAY BACKGROUND 

The intensity of the Cosmic X-Ray Background, the 
integrated measure of the X-ray emission over all red- 
shifts and luminosities, allows us to assess what fraction 
of the population we may be missing from our optically- 
identified point sources in the 2-8 keV band. In addi- 
tion to constraining an undetected population of X-ray 
sources, we can also assess whether our incompleteness 
corrections to account for unidentified sources, are rea- 
sonable. The intensity of the CXRB in the 2-8 keV band 
has been measured with several hard X-ray observat ories 
over the last few decades (e.g.. iMoretti et al]|2003D and 
has been shown to have significant scatter ( ~ 20%) pre- 
sumably due to instrument calibrat ion. Here , wecom- 
pare to the best-fit value from Mor etti et al.l (|2003f ) of 
1.79 X 10"^^ erg cm"^ s~^ deg"^ ^y^^ 2-8 keV band 
that has been converted from the 2-10 keV band based 
on an assumed CXRB spectrum with photon index of 
1.4. 

In Figure 1141 we show the integrated contribution of 
AGN from our best-fit LDDE model (C; solid line) to the 
CXRB as a function of redshift. As shown, we account 
for 52% of the 2-8 keV CXRB at z = 0. We assess the 
contribution of AGN that have X-ray fiuxes above our 
limits though fell out of the sample due to their faint 
optical magnitudes {r',i' > 24) by applying a correction 
based on the assumption that their redshift distribution 
is equivalent to identified sources at similar X-ray fluxes. 
These sources, as shown in the top panels of Figures [3| 
and [4| by the dot-dashed line, mainly occupy the faint 
end of the X-ray flux distribution. Our measure of the re- 
solved fraction of the CXRB now reaches 68% (Fig. [HI 
dashed curve). These results are not surprising since 
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Fig. 14. — Cumulative fraction of the 2-8 keV CXRB as a func- 
tion of redshift. The solid line corresponds to our best-fit LDDE 
model (C) and a corrected version (dashed line) that accounts for 
excluded, faint (r',i' > 24) targets . For comparison, we show 
(d otted curve) the resolv ed fraction for the best fit LDDE model 
of ILa Franca et al.l II2005I) . The vertical line marks half a Hubble 
time. 

IWorslev et aP ^QM) find that ~ 25% of the hard CXRB 
remains to be resolved. This is also consistent with re- 
cent reso lved fration of 79 ± 8% measu red in the 2-8 keV 
band bv iHickox fc Markevitciil (j2006[ ). Both studies at- 
tribute the bulk of this unresolved emission to be from 
absorbed AGN (log Nh ^ 23 - 24 cm^^) ^j^h redshift 
roughly between 0.5 and 1.5. Therefore, at most, we 
may be inaccurately accounting for ~ 10% of the AGN 
population due to our rather simple assu mptions f or the 
opti cally-faint X-ray detections. Both Ue da et al.l (|2003D 
and lLa Franca et al.l (|2005[ ) have measured resolved frac- 
tions close to unity. Other sources of uncertainty may be 
due to the fact that we did not correct our fluxes for X- 
ray absorption and excluded the targets of ea ch ChaMP 
fleld for which a few are bright AGN (K im et al. 2G07b). 
We further illustrate in Figure [TH that the hard CXRB 
is mainly generated by AGN at z < 1 , when the universe 
was more than half its present ag e, as previously de - 
scribed in many works to date (e.g., iBarger et al.ll200l[ ). 

A current limitation of the ChaMP sample is the lack 
of absorbed {Nh > 10^^ cm~^) and luminous {log Lx > 
44) AGN mainly due to the bright o ptical magnitude 
cut fo r optical spectrosc opic foUowup ([Silverman et al.l 
I2005at [Green et al.|[200l . The lack of these AGN may 
not greatly alter our current luminosity function since 
they do not seem to outnumber the unabs orbed (i.e., 
type 1) AGN at these high luminosities. [Ueda et all 
([20031 ) find the ab sorbed fraction at thes e luminosities 
to be ~30% while iLa Franca et al.[ ([2005f ) find an even 
lower fraction of ~ 20%. An accurate assessment is still 
required since the numbers of highly absorbed and lu- 
minous AGN (i.e., type 2 QSOs) are still limited. The 
ChaMP project has initiated a deep optical spectroscopic 
campaign on Gemini to ide ntify the op tically-faint X-ray 
sources (see Fig. 13 of Silve rman et a l. 2005a) that may 
provide a fair sample of these AGN thus removing our 
current selection bias. If we do relax our definition of 
a type II QSO to having a FWHM i 2000 km s"\ as 
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done by many studies (jSzokolv et al.|[2004l : lEckart et al.1 
120061 ). we have found two such QSOs at z > 3 (CX- 
OMP J114036.2+661317, CXOMP J0017 56.6+163006) 
as shown in Figure [71 The SEX SI survey (jEckart et al.l 
12006. 2005l lHarrison et ar2003). with an unprecedented 
sample of 33 luminous AGN lacking broad optical emis- 
sion lines, has demonstrated that deep 8-lOm class spec- 
troscopy is key to identify these optically faint AGN. 

6. COMPARISON WITH OPTICALLY-SELECTED AGN 

It is useful to directly compare our XLF to the latest 
optic al luminosity functions, as also done in recent works 
fe ^lUeda et al.i2003trBarger et al.l2005HRichards et all 
120061) ■ that cover a similar red s hift range. Curr ently, 
the SDSS (iRichards et all [2001 iJiang et all [2006h and 
COMBO-17 (jWolf et al.ll2003f) . that solely detect type 1 
AGN, cover a redshift range out to z ~ 5. We elect to use 
cgs units throughout this work and transform the opti- 
cal luminosity functions to match our XLF. We convert 
the luminosity used for the SDSS in units of absolute 
magnitude {Mf, z = 2) into a monochromatic luminosity 
^K^nnA^ using t he transformation given in Equation 4 of 
IRichards etlHI (|2006[) . The COMBO-17 rest-frame mag- 
nitudes Mi45o are converted to / ° assuming a power- 

2500A 

law slope for the spectral energy distribution {fi, oc v^'^) 
with spectral index (a) equal to 0.5. These monochro- 
matic rest-frame UV luminosities are translated into a 
monochro matic X-ray lumino sity {I2 kev) using Equa- 
tion Ic of ISteffen et al.l ()2006[ ). This relation has been 
shown by many X-ray-to-optical studies of AGN (e.g 
iWilkes et al. [ [l99l iGreen et al.l[l995l : ISteffen et al.l[200l 
to have no redshift depend ence and is applicable out to 
z = 5 (jVignali et al.l[2005h . The slope of th i s rel ation 
is further supported, as shown bv lUeda et al.l (|2003l) . by 
shifting the hard XLF to match the optical at lower red- 
shifts (z < 2.3). We determine the rest-frame 2-8 keV X- 
ray luminosity from /2keV by assuming a pow erlaw SED 
{Le oc E-^) with r = 2.0 ([Reeves fc Turneij [2000). 

We show in Figure [151 a- comparison of our X-ray lu- 
minosity function with the data and analytic model fits 
from the SDSS and COMBO-17 in the four highest red- 
shift intervals mainly to highlight the behavior above and 
below the QSO peak at z 2.5. For the optically- 
selected luminosity functions, we limit the luminosity 
range to that which is covered by these surveys. This 
allows us to illustrate the complementarity of these sur- 
veys and avoid any interpretation based on extrapola- 
tions beyond the dynamic range of each sample. In gen- 
eral, the bright end slope of our XLF agrees with that 
from the optical surveys out to z ~ 5 with the exception 
of the 1.5 < z < 2.0 redshift interval where the SDSS 
model is much steeper even compared to b oth X-rays 
and to COMBO-17 fsee IRichards et al.l [200I . Qualita- 
tively, there does not appear to be any strong evidence 
for a change in the slope of the bright end with redshift 
(see Fig ure [9| for a larger redshift baseline), as s hown 
recentlv ([Richards et al.1l2006HHopkins et al.ll2007i ). At 
lower redshifts (z < 1.5), the slo pe at the bright end i s 
not nearly as steep as measured in lRichards et al.) ([2006D , 
though larger samples of luminous X-ray selected AGN 
are required to justify such a statement. As shown, the 
faint end of the luminosit y function can be probed by 
optically-selected surveys ([Jiang et al.l [20061 : [Hunt et al.l 



l2004f ) of type 1 AGN at z < 3. As expected, their space 
density falls below those found from X-ray selection due 
to the lack of obscured AGN. At z > 3, the luminosities 
of these optically-selected AGN do not fall well below the 
knee. As evident. X-ray surveys can probe moderate- 
luminosity AGN at 3 < z < 4, fainter than the break 
luminosity {log Lx < 44.5), thus allowing a measure 
the slope at the faint end. Some uncertainty remains 
since the contribution of luminous type 2 QSOs is still 
uncertain especially a,t z > 3. However, as previously 
mentioned, recen t constraints on the type 2 QSO popu- 
lation (^25-50% iGilh et all [2001 tend to show hmited 
effect upon our current measure of the XLF even when 
considering th e possible increase in the obscured fraction 
with redshift ([Treister et al.ll2006f) since the luminosity 
dependence is stronger than the redshift dependence. 

7. ACCRETION HISTORY OF SMBHS 

The luminosity function of AGN is one of the key ob- 
servational constraints on the mass buildup of SMBHs 
over cosmological time. With our measure of the XLF 
less biased against obscuration, sensitivity to lower lu- 
minosity AGN and coverage out to high redshift, we 
can reexamine, with very simplistic assumptions, our un- 
derstanding of how the mass accretion rates, Eddington 
masses and the cumulative mass density of SMBHs be- 
haves as a function of redshift. In addition, we can ascer- 
tain whether the well known similarity between the ac- 
cretion rate onto SMBHs and the star formation history 
of galaxies (iBovle &: TerlevichI [1998: Franceschini et aTl 
,1999 ; M erloni et al.ll2004n . from the QSO peak (z ~ 2.5) 
to the present, persists at higher redshifts. 

We determine the distribution of mass accretion rate 
as a function of redshift with the assumption that the 
accretion rate is directly proportional to the bolometric 
luminosity [Lbol] Equation [T^ by a constant mass to 
energy conversion factor trad- 

LbOL — ^radMi^ccC^ (12) 

We derive bolometric luminosities from the rest-frame 
2-8 keV X-ray luminosity using the non-linear relation 
given in Equation 21 of .Marconi et al.. (2004 ) that mainly 
accounts for changes in the overall spectral energy dis- 
tribu tion of AG N as a function of optical luminosity 
(e.g.. ISteffen et al. 2006). We fix the accretion efficiency 
erad to 0.1 fe.g.. lMarconi et al.ll200l . In Figure [TBI we 
plot the accretion rate distribution, using our best-fit 
LDDE model (C), at six redshifts that span the full ex- 
tent of our XLF. Using the same bolometric luminosi- 
ties, we label the top axis with the corresponding Ed- 
dington masses. We stress that these accretion rates 
and masses are highly degenerate. Here, we aim to 
just illustrate the evolutionary behavior of the popula- 
tion in terms of their physical properties rather than 
present well-constrained physical measurements. In light 
of these words of caution, we see that the number den- 
sity of the entire distribution systematically rises with 
cosmic time from z — 4.5 to z = 2.5 (dotted line) that 
may represent the global ignition of mass accretion onto 
SMBHs at early epochs. At redshifts below this peak, 
the number density of AGN with high accretion rates 
(log > 1) or masses i^Edd > 10* M©) essentially 
halts to z ^ 1.5 and then declines rapidly. In contrast. 
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Fig. 15. — Comparison of the X-ray and optical luminosity functions at z > 1.5. The black points an d lines are our data as shown in 
Figu re [9] The op tical luminosity function is overplotted for bot h the SDSS (red line MRichards et al.ll20d3 ), optically-faint SDSS (small red 
dots: IJiang eraU l2006') . and COMBO-17 (dotted red line; Wolf etalj iSoMl . 

the number of AGN with low accretion rates {log M < 1) 
or masses {MEdd ^ 10^ Mq), continues to climb. Ei- 
ther, we are witnessing the emergence of young SMBHs 
(i.e., low mass, high accretion rate) or the remnants of a 
once highly luminous SMBH population (i.e., high mass, 
low accretion rate). Much recent evidence points to a 
luminosity-dependent lifetime fe.g.. Hopkins et "aLll2005l : 
lAdelberger fc Steid cl 2005) that must also be folded in 
to accurately extract physical quantities from a luminos- 
ity function. Evidently, a more careful deconvolution is 
warranted. 

As done in many recent studies, we can measure the 
cumulative mass density of SMBHs as a fu nction of red- 
shift from the integrated light from AGN (|Soltanlll983) 
using the same assumption for the radiative efficiency 
given above. In particular, the local value (z = 0) al- 
lows a direct comparison with complementary techniques 
(e.g., local galaxy luminosity function and a relation be- 
tween black hole mass and its host bulge; integrated light 
from optically-selected QSOs). With the dropoff in the 
number of AGN at z > 3, we do not expect our mea- 
sure of the local value of the mass density to differ from 
recent estimates but aim to illustrate that further con- 
straints on the mean growth of AGN at high redshifts 
are now possible. In Figure [iTl we plot the results us- 
ing our best fit 'mod-PLE' model (D) over an extended 
luminosity (41 < log L2-8 kcV < 47) and redshift range 
(0 < z < 6). We do not account, as done in Section^ for 
the optically-faint X-ray sources not included in our sam- 
ple; their contribution is not trivial to assess since their 




Fig. 16. — Distribution of mass accretion rate onto SMBHs 
at redshifts above and below the peak number density {z = 2.5; 
dashed line) using our best-fit LDDE model (C) and a fixed accre- 
tion efficiency. The top axis shows the Eddington masses assuming 
an equivalent bolometric luminosity. 

redshift distribution is most likely dissimilar to optically- 
bright er AGN at similar X-ray fluxes as demonstrated 
by M artinez-Sansigre et all ([2005). No correction for in- 
trinsic absorption is applied. We have measured the mass 
density of black holes for two extreme values of the radia- 
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tive accretion efficiency tliat is less constrained than the 
bolometric corrections involved. The grey area illustrates 
the range of mass densities with values of the accretion 
efficiency ranging from a Schwarzschild (e^ — 0.06) to a 
Kerr (e^ — 0.3) black hole. The solid line gives the values 
for the widely accepted value e = 0.1 with a mass den- 
sity at z = of 1.64 X 10^ M© Mpc"^. In addition, we 
show the value fo r a constant bolometric correction (40; 
lElvis et al.lll994f ) with e = 0.1. Our estimate of the lo- 
cal mass density in SMBHs is lower than other measures 
(left). Our value is 48% of the the local value, deter- 
mined from the galaxy luminosity functi on and velocity 
disperson, of 3.4^0-5 x 10^ Mq Mpc'^ (jMarconi et al.l 
[2OOI and is in closer agreement to the va lues obtained 
from previous hard X-ray-selected A GN (B arger et al.l 
120051 iLajranca et aLllgoalUeda et~al. 2003) . In partic- 
ular, iBarger et al.l ( 20051 ) find a similar local black hole 
mass density (2.1 x 10^ Mq Mpc""^) when using a bolo- 
metric correction fa ctor (35) close to that employed here. 
I Barger et a l. ( 2005) point out th at the mass density mea- 
sured bv lYu fc Trcmaini (|2002D is possibly too high due 
to an extrapolation of the optical luminosity function at 
the faint end that has now been shown to be shallower 
(jCroom et al.ll2004D . The fact that we are not accounting 
for optically-faint X-ray sources, resolving ^70% of the 
hard CXRB, and do not account for Compton-thick AGN 
could contribute to the discrepancy with mass density es- 
timates using local galaxies. As is evident in Figure [T7] 
right, most (81%) of the growth of SMBHs occured by 
the time the universe was half (z ^ 0.7) its present age 
with minimal contribution at z > 3. 

7.1. Co-evolution of SMBHs and star-forming galaxies 

The order of magnitude decline in the co-moving mean 
emissivity of AGN and star formation history of galax- 
ies from z ~ 1.5 to the present has initiated many early 
postulates for a co-evolution s cheme of SMBHs and the 
galaxies in which they reside. iBovle fc TerlevichI ()1998l ) 
demonstrated that the optical QSO luminosity function 
exhibited a similar beh avior to the sta r formation rate of 
field galaxies. Franccs chini et al.l (fl99 9). with an X-ray 
selected sample of AGN, further investigated this cor- 
relation by demonstrating that the more luminous QSO 
evolve closely with the most massive galaxies (i.e., E/SO) 
while the lower-luminosity AGN track the star-forming, 
field galaxy population. At the time, these studies were 
mainly limited to 2 < 3 due to the availability of mea- 
sured SFRs of galaxies. 

Armed with our XLF of AGN out to z '--^ 5 and 
the latest SFRs of high redshift galaxies, we can revisit 
this correlation. In the last few years, new samples of 
high redshift galaxies has enabled the ev olution of the 
SFR to be extended out t o z ~ 6 fe.g.. iBunker et al.l 
120041: iGiavalisco et al.ll2"00l iBouwens et al.ll2007D . Most 
studies have selected these high redshift galaxies by Ly- 
man break techniques (e.g.,^^Stcidel et al. 1999) or narrow 
emission line searches fe.g.. IRhoads et al.ll2003D . With 
the advent of wide area, multi-slit spectrographs on 8- 
10m class telescopes, flux-limited surveys are able to 
probe a sizeable volume and identify significant numbers 
of high redshift galaxies without preset color selectio n 
criteria or emission line strengths (|Le Fevre et al.l[2005l ). 

We have measured the mass accretion rate per co- 
moving volume (units of Mq yr~^ Mpc~^) as a function 
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Fig. 17. — Integrated mass density of SMBHs. Left The dark 
grey region shows our cumulative values using our best fit mod- 
PLE model (D), a luminosity-dependent bolometric correction, and 
a range of accretion efficiencies of 0.06 < Cr < 0.3 with a solid line 
marking the Cr = 0.1 case. We plot with a dashed line the results 
using a constant bolometric correction and er = 0.1. For com- 
parison, we have highlighted the 2 = mass densities (horizontal, 
dash ed lines) measured using the optical QSO luminosity func- 
tion llYu fc Tr emainc 2002) and hard XLF (Marconi et al. 2004) of 
lUeda et ain 2003). The dotted line marks the mass density as a 
function of redshift using the hard XLF of La Franca et al. (2005). 
Right We show the cumulative mass density, identical to the case 
on the left marked by a solid line, as a function of lookback time 
normalized by our local value. 



of redshift using our analytic models (C, D) for all AGN 
with L2-S kcV > 10^^ erg s^^. Bolometric corrections 
and conversion factors to accretion rates are used as de- 
scribed in the previous section. In Figure 1181 we plot 
the results compared to the dust-corre c ted S FRs (black 
data points) as compiled by iHopkinsI (l2004f). We add 
the la test measurements at z > 4 from IBouwens et al.l 
()2007t ) as shown by the open circles. The mass accretion 
rates have been scaled up by a factor of 5000. Up to 
z ~ 2, there is an overall similarity as previously elab- 
orate d in many s t udies to date. The discrepancy noted 
by Ha iman et all ()2004[ ) at z < 2 appears to be removed 
with the inclusion of obscured AGN from X-ray selected 
surveys. The peak in the black hole accretion rate den- 
sity now occurs at z ~ 1.5 rather than z ~ 2.5 as mea- 
sured with optically-selected QSO samples. Also noted 
in many investigations, there is a divergence at z > 2 
between the two growth rates, with a faster decline of 
the A GN population even w hen comparing to the recent 
SFRs (jBouwens et al.l [20071) that show a significant de- 
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Fig. 18. — Comparison of the star formation history of galaxies 
(data points) to the mass accretion rate ( Mq yr~^ Mpc~^) of 
SMBHs using the Mod-PLE (solid line) and LDDE (dotted line) 
models. SFRs from the compilation of iHopkina (2004) are marked 
by the small filled dots wh ile the open circles are recent measure- 
ments at 2 > 4 (Bouwens et al.|[2007I V The mass accretion rates 
have been scaled up by a factor of 5000. The divergence at 2: > 4 of 
the two curves representing the accretion history of SMBHs high- 
lights the poor constraints on the X-ray luminosity function at 
early epochs. 

cline between 4 < z < 6. Here wc find a decrease in the 
accretion rate onto SMBHs by a factor of 3-4 compared 
to the SFRs at z ~ 4.5. This may represent a delay be- 
tween the formation of stars and fully matured SMBHs 
in the early universe possibly due to either the cooling 
times required for gas to be avai lable for accretion o r an 
insufhcient merger history (See iHaiman et al.|[200^ . for 
further discussion) . 

8. DETECTION OF HIGH REDSHIFT AGN IN CHANDRA 
SURVEYS 

Even though we have extended the 2-8 keV luminos- 
ity function beyond z ~ 3, there remains much uncer- 
tainty about the behavior of the luminosity function be- 
low the break in the 3 < z < 4 redshift range and at 
all luminosities at higher redshifts. Since Chandra has 
now completed or will undertake some new surveys that 
cover a wider area and at moderate depths, we can make 
some predictions of their return of high redshift AGN. We 
have computed the numbers of AGN in two redshift inter- 
vals and two luminosity ranges that basically span below 
(42 < log Lx < 44) and above (44 < log Lx < 46) the 
break luminosity. We give these values for various sur- 
veys using our best-fit LDDE and mod-PLE (in paren- 
thesis) models in Table [S] The estimates between the 
two models are similar at high luminosities but differ at 
the low luminosity end, which becomes even more pro- 
nounced at the highest redshifts (see Figure [H]). 

We see that the ChaMP will provide a total of 
46 AGN at z > 3 based on an area coverage of 
2 square degrees that is similar to the numbers ex- 
pected from the latest wider area surveys (i.e., EGS, 
Chandra/ COSMOS) . As mentioned earlier, these surveys 
are complementary to ChaMP since they reach greater 
depths over a smaller area to detect AGN at lower lumi- 
nosities as exemplified here. We note that the mod-PLE 



model is most likely overpredicting the numbers of lower 
luminosity {log Lx < 44) AGN based on the observed 
number in the CDF-N and CDF-S. It total, we expect a 
final sample of over ~195 AGN at z > 3 with 63 above a 
redshift of 4. Better constraints on the faint end of the 
luminosity function in the 3 < z < 4 redshift interval can 
be achieved by either the E-CDF-S, EGS, or cCOSMOS 
with a sample of ~ 15 — 25 sources per survey though an 
accurate measure of the slope will require the combined 
sample that is especially true at higher redshifts. An ac- 
curate assessment of the luminosity function at z > 4 
will definitely require a merger of all available catalogs. 

9. SUMMARY 

We have presented an extension of the hard (2-8 keV) 
X-ray luminosity function of AGN up to z ~ 5. The 
ChaMP effectively covers a wide area (1.8 deg^) at suffi- 
cient depths (/o. 5-2.0 kov ^ 10~^^ erg cm~^ s~^) to sig- 
nificantly improve the statistics of luminous (Lx > 10*** 
erg s^^) AGN at high redshift. In total, we have amassed 
a sample of 682 AGN with 31 at z > 3. The addition of 
lower luminosity AGN from the Chandra Deep Fields is 
instrumental to characterize the faint end slope at z > 1 . 
We have corrected for incompleteness (i.e., fraction of 
sources without a redshift) as both a function of X-ray 
flux and optical magnitude, as dictated by the limita- 
tions of optical spectroscopic foUowup. Significant opti- 
cal foUowup is still required to accurately account for the 
obscured population, especially at high redshifts. 

We have measured the hard XLF with both a binned 
(l/Va) and an unbinned method that fits an analytic 
model to our data using a maximum likelihood tech- 
nique. We found that the luminosity function is simi- 
lar to that foun d in p revious studie s (lUeda et ani2003t 
iLa Franca et al.l 120051 : fBarger et al.1l2005D up to^ = 3, 
with an evolution dependent upon luminosity. At higher 
redshifts, there is a significant decline in the numbers of 
AGN with an evolution rate similar to that found with 
optically-selected QSO samples. We further show that 
the strong evolution above the cutoff redshift may cause 
the LDDE model to underpredict the number of low lu- 
minosity AGN at z > 2. A PLE model with a faint-end 
slope dependent on redshift agrees better with the binned 
(1/Va) data at z ~ 2.5 though it may overpredict the 
number of faint AGN at higher redshifts. We highlight 
the need to improve the statistics of both high redshift 
AGN at z > 4 and lower luminosity (i.e., below the break 
in the luminosity function) AGN at z > 3. Such improve- 
ments are feasible from our predictions of the numbers 
AGN that will be found in the latest Chandra surveys 
(E-CDF-S, EGS and COSMOS) 

Our new luminosity function accounts for 52% of 
the 2-8 keV Cosmic X-ray Background. The integrated 
emission from these AGN give a z = mass density 
of SMBH of 1.64 X 10^ Mg Mpc'^, lower than other 
published values using X-ray selected AGN samples and 
the local value measured using a galaxy luminosity func- 
tion and a bulge-velocity relation, possibly due to unac- 
counted optically-faint AGN {r',i' > 24), and Compton- 
thick AGN. Further, a comparison of the mean mass ac- 
cretion rate of SMBHs to the star formation history of 
galaxies out to z ~ 5 shows the familiar co-evolution 
scheme up to z 2 and a divergence at higher redshifts 
with perhaps star formation preceeding the formation of 
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SMBHs. 
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TABLE 1 

AGN SELECTED FROM VARIOUS X-RAY SURVEYS 



Name 


Telescope 


Hard band selection (z 


<3) 


Soft band selection ( 


^ > 3) 






Depth 


r' 


Number 


Depth 


i' 


Number 






(erg cm~^ 


(mag) 


of AGN 


(cgs) 


(mag) 


of AGN 


ChaMP 


Chandra 


2.7 X 10-15 


22.0 


273 


1.0 X IQ-i"' 


22.0 


13 


CLASXS 


Chandra 


2.3 X 10-15 


24.0 


103 


5.6 X 10-1*5 


24.0 


3 


CDF-N 


Chandra 


6.3 X 10-18 


24.0 


104 


1.0 X 10-16 


24.0 


8 


CDF-S 


Chandra 


6.3 X 10-16 


24.0 


95 


1.0 X 10-16 


24.0 


2 


Lockman Hole 


XMM-Newton 








1.0 X 10-15 


24.0 


5 


AMSSn 


ASCA 


2,52 X 10-1'^ 


22.0 


76 









TABLE 2 
ChaMP fields 



Obs. ID 


ACIS 


PI Target 


Exposure 
time (ksec) 


520 


I 


MS0015.9-I-1609 


61.0 


913 


I 


CLJ0152.7-1357 


34.6 


796 


I 


SBS0335-052 


47.0 


624 


S 


LP944-20 


40.9 


902 


I 


MS0451. 6-0305 


41.5 


914 


I 


CLJ0542.8-4100 


48.7 


377 


s 


B2 0738+313 


26.9 


2130 


s 


3C207 


30.0 


419 


s 


RXJ0911. 4+0551 


30.0 


839 


s 


3C220.1 


20.0 


512 


s 


EMSS1054.5-0321 


75.6 


363 


s 


PG1115+080 


30.0 


536 


I 


MS1137.5+6625 


114.6 


874 


I 


1156+295 


75.0 


809 


s 


Mrk237X 


50.0 


541 


I 


V1416+4446 


29.8 


800 


s 


CB 58 


50.0 


546 


I 


MS1621. 5+2640 


30.0 


830 


s 


Jet of 3C390.3 


23.6 


551 


I 


MS2053.7-0449 


42.3 


928 


s 


MS2137.3-2353 


29.1 


431 


s 


Einstein Cross 


21.9 


918 


I 


CLJ2302.8+0844 


106.1 


861 


s 


Q2345+007 


65.0 



TABLE 3 
ChaMP statistics 



X-ray 
band 


X-ray flux 
limit* 


Optical 
limit 


All 


Numbers 
Spectra 


Ids'" 


Id fraction 
All observed 


Hard 


2.7 X 10-15 


none 


793 


456 


328 


41% 


72% 






r' < 24.0 


598 


444 


325 


54% 


73% 






r' < 22.0 


390 


332 


289 


74% 


87% 


Soft 


4.0 X 10-16 


none 


1125 


624 


435 


39% 


70% 






i' < 24.0 


834 


605 


430 


52% 


71% 






i' < 22.0 


512 


429 


373 


73% 


87% 



^ units of erg cm ^ s ^ 

^ Quality flag equal to two or three. 



TABLE 4 

Hard XLF Best-Fit Model Parameters 



Model 


Type 


log Ao 


7I 


72 


log Lo 


cl 


e2 


Zc 


log La 


a 


Redshift range 


A 
B 
C 
D 


PLE 
LDDE 
LDDE 
Mod-PLE 


-5.332±0.015 
-6.077±0.015 
-6.163±0.015 
-5.238±0.015 


2.80±0.20 

9 it;+0.42 
'^■'-•^-0.12 

2.15'" 

2 76+0-21 
^•'"-0.19 


0.77±0.05 
L10±0.13 
LIO'" 
0.42±0.05 


44.92±0.10 
44.33±0.10 
44.33''''^ 
44.88i°:^^ 


1.69±0.12 
4.00±0.28 
4.22t°-i? 


-0.87±0.6 

-1.5=^ 

r,7 + 0.31 

-0.34 

-8.18±0.55 


1.9=" 
1.9=" 

2.0{Fixod) 


44.6=' 
44.6=^ 


n 017+0.020 
U.Ol ( _o.027 

0.333±0.013 

1 Q4+0.11 


0.2 < z < 3.0 
0.2 < z < 3.0 
0.2 < 2 < 5.5 
0.2 < 2 < 5.5 



^ Fixed parameters to those best fit values from fTjeda et all d2003f) 
^ Fixed to the best fit value in Model B. 
^ For this case, L* — Lo 



TABLE 5 

Predicted number of high redshift AGN in complete surveys 



Survey 


3 < 


2; < 4 


4: < z 


< 6 


Total 




log Lx < 44 


log Lx > 44 


log Lx < 44 


log Lx > 44 




ChaMP 


3(6) 


31 (41) 


0(0) 


12 (6) 


46 (53) 


CDF-N 


4 (16) 


3 (5) 


2(9) 


2 (2) 


11 (32) 


CDF-S 


4 (12) 


3 (4) 


2(6) 


2 (2) 


11 (24) 


E-CDF-S 


11 (35) 


8 (13) 


4 (18) 


6 (5) 


29 (71) 


EGS 


12 (38) 


16 (25) 


4 (13) 


11 (8) 


43 (84) 


cCOSMOS 


15 (45) 


22 (35) 


4(14) 


14 (12) 


55 (106) 


Totals 


49 (152) 


73 (114) 


16 (60) 


47 (35) 


195 (370) 







CD 

S- 

CD 



to 
cn 



